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In previous studies the pressure dependence of the magnetic ordering temperature To of Dy was 
found to exhibit a sharp increase above its volume collapse pressure of 73 GPa, appearing to reach 
temperatures well above ambient at 157 GPa. In a search for a second such lanthanide, electrical 
resistivity measurements were carried out on neighboring Tb to 141 GPa over the temperature range 
3.8 - 295 K. Below Tb’s volume collapse pressure of 53 GPa, the pressure dependence To{P) mirrors 
that of both Dy and Gd. However, at higher pressures To{P) for Tb becomes highly anomalous. 

This result, together with the very strong suppression of superconductivity by dilute Tb ions in 
Y, suggests that extreme pressure transports Tb into an unconventional magnetic state with an 
anomalously high magnetic ordering temperature. 


I. INTRODUCTION 

The magnetic ordering temperatures To of Gd and 
Dy have been recently shown to track each other in a 
highly non-monotonic fashion as a function of pressure 
to ^ 70 GPa; at higher pressures they deviate markedly. 
To for Dy rising rapidly to temperatures well above am¬ 
bient at 157 GPa.Ii^ Parallel experiments on dilute mag¬ 
netic alloys of Gd and Dy with superconducting Y sug¬ 
gest that for pressures above ^ 70 GPa Dy is trans¬ 
formed from a magnetically conventional lanthanide into 
one with an unconventional magnetic state with marked 
deviations from de Gennes scaling,a state perhaps gov¬ 
erned by Kondo physics, indicating that the Dy ion is 
nearing a magnetic instability.® An alternate explana¬ 
tion is that the strong enhancement of To in Dy arises 
through changes in the crystalline electric field at ex¬ 
treme pressure.® Analogous studies on additional lan¬ 
thanides are recommended to help identify the origin of 
this anomalous behavior. 

The lanthanide Tb, which lies between Gd and Dy in 
the periodic table, has one fewer 4/ electron than Dy, 
and is probably less stable magnetically than Dy due 
to its direct proximity to Gd, by far the most stable 
of all magnetic lanthanides. Tb orders antiferromagnet- 
ically (AFM) at To ~ 230 K followed by a ferromag¬ 
netic (FM) transition at To ~ 220 K.® Both transition 
temperatures initially decrease rapidly with pressure at 
the rate -10 to -12 K/GPa, but above ^ 7 GPa nei¬ 
ther transition can be clearly detected in either the ac 
or dc magnetic susceptibility.®® The disappearance of 
the ordered moment in the susceptibility measurement 
indicates a transition to either an AFM or paramag¬ 
netic state above 7 GPa. Electrical resistivity studies 
should reveal which scenario is correct since both FM 
and AFM order normally lead to a distinct kink in the 
temperature dependence of the resistivity. However, re¬ 
cent resistivity and neutron diffraction experiments on 
Tb find that the FM transition decreases with pressure 
at the rate -16.7 K/GPa to 3.6 GPa® that the transi¬ 


tion could no longer be resolved above 3.6 GPa may be 
due to appreciable pressure-gradient broadening in the 
cell which contained no pressure medium. That mag¬ 
netic order in Tb disappears above 7 GPa seems highly 
unlikely since both x-ray absorption near-edge structure 
(XANES) and non-resonant x-ray emission spectroscopy 
(XES) measurements detect no chaise in Tb’s valence 
to 65 GPa and 70 GPa, respectively.® In fact, the XES 
studies show that Tb retains its strong, highly localized 
magnetic moment (J = 6) to at least 70 GPa.® 

In this paper we present the results of dc electrical re¬ 
sistivity measurements on Tb over the temperature range 
3.8 - 295 K to pressures as high as 141 GPa, well above 
the pressure of 53 GPa where Tb suffers a 5% volume 
collapse at the phase transition from hexagonal hR24 to 
body-centered monoclinic (bcm).® Magnetic order is in¬ 
deed observed in Tb for pressures above 7 GPa. In fact, 
to 53 GPa To(T) follows nearly the same highly non¬ 
monotonic pressure dependence found earlier in Gd and 
Dy,® but deviates markedly at higher pressures. As the 
applied pressure passes through 53 GPa, To(P) for Tb 
first decreases, but then begins to increase rapidly above 
80 GPa. As suggested for Dy,® extreme pressure appears 
to transport Tb into an unconventional magnetic state 
with an anomalously high magnetic ordering tempera¬ 
ture, well above that anticipated from conventional de 
Gennes scaling. 


II. EXPERIMENTAL TECHNIQUES 

Resistivity samples were cut from a Tb ingot 
(99.9% pure, Material Preparation Genter of the Ames 
Laboratorj®!) Xo generate pressures well beyond the 
volume collapse pressure of Tb at 53 GPa, a dia¬ 
mond anvil cell (DAG) made of CuBe alloy was used.l^ 
Two separate high-pressure experiments were carried out 
where pressure was generated by two opposing diamond 
anvils (1/6-carat, type la) with 0.35 mm diameter culets 
beveled at 7° to 0.18 mm central flats. 
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The Re gasket (6-7 mm diameter, 250 fim thick) 
was preindented to 30 /rm and a 80 /im diameter hole 
electro-spark drilled through the center. The center sec¬ 
tion of the preindented gasket surface was hlled with a 4:1 
cBN-epoxy mixture to insulate the gasket and serve as 
pressure medium. The thin Tb sample was then placed 
on top of four thin (4 /rm) Pt leads for a four-point dc 
electrical resistivity measurement. In an attempt to min¬ 
imize the effect of the pressure gradient across the sample 
in this non-hydrostatic pressure environment, in run 1 an 
elongated sample (dimensions ^8x80x3 ^m^) was used 
with the two voltage leads spaced only 5 /rm apart (see 
inset to Fig. 1(a)). In run 2 all four Pt leads were placed 
near the corners of the square-shaped sample (dimen¬ 
sions ~30x30x5 fj,m) (see inset to Fig. 1(b)), as in the 
previous resistivity measurements on Dy.^^ However, from 
the temperature-dependent resistivity data the pressure 
gradient was estimated to be approximately the same in 
both runs. Further details of the non-hydrostatic high 
pressure resistivity technique are given in a paper by 
Shimizu et aim 

A He-gas driven membrane was utilized to change pres¬ 
sure at any temperature.^ The value of t he p ressure 
was determined using both the fluorescenceHl from a 
small ruby sphere positioned at the edge of the sample 
and the frequency shift of the diamond vibron via Ra¬ 
man spectroscopy.^ The ruby pressure was determined 
at both ambient temperature and a temperature within 
20 K of To; the vibron pressure was determined only at 
ambient temperature. The values of the pressure given 
are averaged over the sample to an estimated accuracy 
of ±10%. In these experiments temperatures from 3.8 
K to 295 K were available using an Oxford flow cryo¬ 
stat. All measurements shown in this paper were car¬ 
ried out with increasing pressure; diamond anvil failure 
at the highest pressure ended the experiment. Further 
experime ntal detai ls of the DAC and cryostat are given 
elsewhere! ^ * ^^ * ^^ * ^^ 


III. RESULTS OF EXPERIMENT 

The present resistivity studies on Tb were carried out 
in two separate experiments. In Fig. 1(a) the electrical 
resistance R{T) from run 1 is plotted versus temperature 
at 18 different pressures to 135 GPa. The results from run 
2 are shown in Fig. 1(b) and span the pressure range 2 - 
141 GPa with 17 values. The onset of magnetic ordering 
is identified by the kink in the R{T) dependence clearly 
seen near 200 K at 2 GPa, the lowest pressure in each 
run. The kink in R{T) upon cooling marks the beginning 
of the suppression of spin-disorder scattering Rsd{T) as 
magnetic ordering sets in.^ At higher pressures this kink 
broadens somewhat into a ’’knee” due to an increasing 
pressure gradient across the sample, but remains clearly 
visible to ^115 GPa. 

In Fig. 2 selected data from Fig. 1(b) are replotted 
but shifted vertically for clarity so that no curves inter¬ 


sect. In this graph the red line through the data above 
the knee gives the phonon contribution RphiT) to the 
total measured resistance R{T) estimated in the same 
manner as in our previous work on Dy,!^] as outlined in 
the next paragraph. The paramagnetic state of Tb yields 
the relatively flat region of R{T) at higher temperature 
where the red (phonon) line overlaps the data. Where 
the red line begins to separate from the data marks the 
onset of magnetic ordering in some region of the sample. 
Because of the pressure gradient across the sample, other 
regions of the sample will have a lower onset temperature, 
thus broadening the kink into a knee. The intersection 
of the phonon resistance (red curve) with the red low- 
temperature tangent curve defines the temperature Tx in 
Fig. 2. 

The total measured resistance is the sum of three 
terms, RiT) = Rd-\- RphiT) ± Rsd{T), where Rd is the 
temperature-independent defect resistance, Rph{T) the 
temperature-dependent phonon resistance, and RsdiT) 
the temperature-dependent spin-disorder resistance. At 
temperatures where there is no magnetic ordering in the 
sample, RgdiT) is independent of temperature. Above 
the onset temperature of the knee, the temperature de¬ 
pendence of R{T) is, therefore, due solely to that of 
RphiT). The temperature dependence of the phonon re¬ 
sistance is visible over the widest temperature range at 
that pressure (16 GPa in Fig. 2) where the knee begins 
at the lowest temperature. We extrapolate this depen¬ 
dence to 0 K in the temperature region below the knee 
to yield the temperature-dependent function (T), the 
estimated phonon resistance at 16 GPa in run 2. In run 1 
the data at 18 GPa were used in the same way to obtain 
Rp®(T). Since the functional dependence of Rph{T) on 
temperatures above Tx is seen in Fig. 2 to change only 
slowly with pressure, we estimate RphiT) for the other 
pressures in run 2 by simply multiplying the function 
RphiT) by a ’’phonon factor” a chosen such that for tem¬ 
peratures above the knee the quantity i?(T) — aR^iT) 
becomes temperature independent for T > T^. The val¬ 
ues of a required are listed in Table 1 at all pressures in 
run 2 to 141 GPa. For pressures of 120 GPa and above, 
the knee in RiT) apparently begins above 295 K, so that 
a can no longer be estimated directly from the resistance 
data. For P > 120 GPa, therefore, the value a = 0.41 
is assumed in Table 1 for run 2 and a = 0.69 in run 1. 
This is admittedly an oversimplified way to estimate the 
phonon contribution, but is superior to the assumption 
made in an earlier study that for many lanthanides Rph 
is simply a linear function of temperature to 0 

In Fig. 3 the extracted spin-disorder resistance 
RsdiT) = RiT) — aRphiT) — Rd is plotted for pressures 
5, 35, and 86 GPa in run 2. The saturation (maximum) 
value of the so obtained spin-disorder resistance in 

the paramagnetic phase at each pressure is listed in Table 
1 for run 2. At 86 GPa, for example, — 152 mU as 
seen in Fig. 3. A similar procedure was used to obtain 
in run 1. 

As illustrated in Fig. 3, the average magnetic order- 
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FIG. 1. (color online) Resistance of Tb versus temperature to 295 K for (a) ruu 1 and (b) run 2 at various pressures. Insets 
show photograph of (a) elongated sample in run 1 and (b) square sample in run 2. 


ing temperature Tq in the Tb sample is estimated from 
the point of intersection of two straight (red) lines, a 
horizontal line for temperatures above the onset of the 
knee, and a line tangent to Rsd{T) at lower tempera¬ 
tures. For 86 GPa it is seen that Tq ~ 69 K. This 
temperature differs by only 1 K from Tx ~ 68 K, the 
intersection point of the phonon resistance and the low- 
temperature tangent lines in Fig. 2 at the same pressure. 
Here we regard Tq to be the average magnetic ordering 
temperature (in our previous paper on Dy, Tx was used as 
the ordering temperatur^. Since the pressure gradient 
leads to a variation in the value of the magnetic order¬ 


ing temperature across the sample, we define the “maxi¬ 
mum” ordering temperature as the temperature at 

which the spin-disorder resistance has decreased by 1%. 
In Fig. 3 it is seen that ~ 113 K at 86 GPa. If 
dTo/dP > 0, gives the value of the magnetic order¬ 

ing temperature at the center of the cell (sample) where 
the pressure is highest. In Fig. 3 it is seen that 
lies 44 K higher than To at 86 GPa. All values of To and 
jimax 2 are listed in Table 1. 

In Fig. 4 To and T™“ are plotted versus pressure to 
141 GPa for runs 1 and 2 on Tb; values for P > 120 GPa 
are estimated using a procedure from Ref. [TJ as outlined 
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FIG. 2. (color online) Selection of resistance versus tempera¬ 
ture curves for Tb from run 2 in Fig. 1(b) where, except at 2 
GPa, the curves have been shifted vertically for clarity. Red 
lines with small positive slope give temperature dependence 
of phonon resistance for T > Tx except at 16 GPa where the 
phonon resistance extends to 0 K (see text). 


below. Where they can be compared, the present results 
are in reasonable agreement with earlier ac magnetic sus¬ 
ceptibility measurements of Jackson et al. to 6.3 GPa.^ 
The pressure dependence To{P) at higher pressures is 
seen to be highly non-monotonic, presumably in response 
to multiple structural phase transitions^ (see top of the 
graph). Note that the phase boundaries were determined 
from x-ray diffraction studies at ambient temperature 
and may shift somewhat as the temperature is lowered. 

A comparison of T'o(P) for Tb from Fig. 4 to compa¬ 
rable graphs for Gd and Dy in Ref. [T] reveals a remark¬ 
able similarity to 53 GPa, the pressure at which the 5% 
volume collapse in Tb occurs.*^ Also plotted in Fig. 4 
are the values of for Tb given by the upper error 

bars connected to the values of To at each pressure by 



FIG. 3. (color online) Spin-disorder resistance Rsd{T) ver¬ 
sus temperature at three pressures from run 2. The phonon 
Rph(T) and defect Rd resistances have been subtracted off. 
The average magnetic ordering temperature To is defined by 
intersection point of two tangent lines. gives temper¬ 

ature at which spin-disorder resistance has decreased by 1% 
(see text). 


a light (blue) vertical line. Particularly intriguing is the 
decrease in Tq following the hR24 to body-centered mon¬ 
oclinic (bcm) transition at 53 GPa,l^ followed by a rapid 
increase above 80 GPa. In contrast to the findings for 
P < 53 GPa, at higher pressures To(P) for Tb thus dif¬ 
fers significantly from that found earlier for either Gd or 
Dy.Ii] Plotted versus relative volume VjVo, the increase 
of To above 80 GPa for Tb is found to be much more 
rapid than the initial decrease of To to 6.3 GPa. A sim¬ 
ilar result was found for Dy.® Extrapolating To versus 
V/Vo for Tb linearly to V/Vo = 0.40 (141 GPa), yields 
the values To « 250 K and T/^^^ « 350 K. 

We now attempt a more quantitative estimate of the 
pressure dependence of To, T™“, and in the pres¬ 
sure range above 116 GPa where the onset of the knee 
appears to lie at or above ambient temperature. We first 
consider the spin-disorder resistance Rsd{T) at pressures 
P < 120 GPa. The first step is to normalize Rsd{T) to its 
value at 295 K, yielding the relative spin-disorder resis¬ 
tance Psd(T)/P^“ plotted versus logT for data at 106, 
112, and 116 GPa in Fig. 5. Since at the higher pressures 
of 120, 126, and 141 GPa the onset of magnetic order¬ 
ing appears to lie above the temperature range of the 
present experiments (295 K), one cannot determine the 
value of Rsd in the paramagnetic phase, nor To or T™“, 
directly from the resistance data. However, noticing that 
over much of the temperature range the Rsd/T) curves 
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FIG. 5. (color online) Relative spin-disorder resistance 
Rsd{T) /versus logT. For data at 120, 126, and 141 GPa, 
RTd^ is estimated by adjusting slope of temperature depen¬ 
dence to match that at 106, 112 and 116 GPa (see text). 
From relative horizontal shifts of the curves the pressure- 
dependence of the magnetic ordering temperature To is es¬ 
timated for 120, 126, and 141 GPa (see text and Table 1). 


FIG. 4. (color online) Average magnetic ordering tempera¬ 
tures To of Tb versus pressure: (-I-) earlier susceptibility stud¬ 
ies to 6.3 GPa with slope dTo/dP — —11 K/GPaPpresent re¬ 
sistance measurements to 141 GPa from (a) run 1, (•) run 2. 
“error bar” connected vertically to each value of To gives maxi¬ 
mum ordering temperature at that pressure. Open sym¬ 
bols indicate extrapolated values (see text). Vertical dashed 
line marks pressure of volume collapse for Tb at 53 GPa. 
Grystal structures for Tb are given at top of graph.^^ R^d^ 
versus pressure is plotted in lower part of figure from run 2 
where it is seen to roughly track To(P). 


age magnetic ordering temperature To has increased from 
148 K to 203 K and the maximum ordering temperature 
jmax 249 K to 341 K, values close to those obtained 
in the linear extrapolation above. 

In our previous work on Dy the spin-disorder resistance 
in the paramagnetic state was found to approxi¬ 

mately track the magnetic ordering temperature To as a 
function of pressure. This same result is seen in Fig. 4 
to hold for Tb. 


for 106, 112, and 116 GPa are approximately parallel on 
the logT plot in Fig. 5, we divide the Rsd{T) data for 
P > 120 GPa by that factor which results in curves par¬ 
allel to those at the lower pressures, as seen in Fig. 5. 
We identify this factor as the value of the temperature- 
independent spin-disorder resistance in the param¬ 

agnetic phase, as listed in Table 1. This procedure is tan¬ 
tamount to assuming that Rsd = Rsd{T/To) for P > 106 
GPa. The change in the value of the magnetic ordering 
temperatures Tq and can now be estimated from 

the shift of the Rsd{T)/R^^^ curves along the logT axis. 
The resulting values of To(T) and T™®'^(T) are given in 
Table 1 and in Fig. 4 for all pressures in runs 1 and 
2 as the open triangles and circles, respectively. From 
this analysis we infer that from 116 to 141 GPa the aver- 


IV. DISCUSSION 

We now seek to identify the mechanism(s) responsible 
for the highly non-monotonic dependence of Tb’s mag¬ 
netic ordering temperature To on pressure. First we fo¬ 
cus on the pressure region below 53 GPa, the pressure at 
which Tb suffers a 5% volume collapse. Since the pres¬ 
sure dependence of To is so similar for Tb, Gd and Dy in 
this pressure range, a common mechanism seems likely. 

For a conventional lanthanide metal with a stable 
magnetic moment, the magnetic ordering temperature 
To is expected to scale with the de Gennes factor 
{g— 1)^ Jt(Jt -I-1), modulated by the prefactor J'^N{Ef), 
where J is the exchange interaction between the 4/ ion 
and the conduction electrons, N{Ep) the density of states 
at the Fermi energy, g the Lande-g factor, and Jt the to- 
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tal angular momentum quantum numberP Since the de 
Gennes factor is constant under pressure, unless the mag¬ 
netic state becomes unstable and/or a valence transition 
occurs, the marked similarity between the highly non¬ 
monotonic pressure dependences of To for Tb, Dy and Gd 
to 53 GPa likely originates in the pressure dependence of 
the prefactor J‘^N{Ep), facilitated by a series of nearly 
ide ntical structural phase transitions in Tb,l^ Dyp^ and 
Qd) 2 i | 22 | "pjjggg phase transitions are likely driven by in¬ 
creasing 5c?-electron occupation with pressure.^ Indeed, 
electronic structure calculations for Dy suggest that its 
large negative initial pressure derivative dTo/dP ^ —6. 7 
K/GPa results from a strong decrease in 

We now consider the pressure region P > 53 GPa 
where the pressure dependence To (P) for Tb is highly 
anomalous, deviating markedly from that of the model 
conventional lanthanide Gd to at least 127 GPaP The 
absence of magnetic instabilities in Gd, even at extreme 
pressures, is expected since the local magnetic state of Gd 
with its half-filled 4/^ shell is the most stable of all ele¬ 
ments, its 4/^ level lying ~ 9 eV below the Fermi levelf^ 
Why is To{P) in Tb anomalous for P > 53 GPa? A long¬ 
standing strategjEZEH to probe the magnetic state of a 
given ion is to alloy this ion in dilute concentration with 
a host superconductor and determine APc, the degree of 
suppression of the host’s superconducting transition tem¬ 
perature. Yttrium (Y) is the ideal host superconductor 
for Tb since the character of its spd-electron conduction 
band closely matches that of the heavy lanthanides, Y 
even exhibiting nearly the same sequence of structural 
transitions under pressure.!^ One may thus anticipate 
that changes in the magnetic state of the Tb ion in the 
dilute alloy will be mirrored in the changes occurring in 
the magnetic state of Tb metal. 

The efficacy of this strategy is supported by previ¬ 
ous studies of the pressure dependences To{P) for Dy 
metal and ATc{P) for Y(Dy) where both experience a 
dramatic enhancement beginning just above the pressure 
of Dy’s volume collapse at 73 GPa.® It was argued that 
this anomalous behavior might be the result of the Dy 
ion exhibiting Kondo physics at elevated pressures where 
both To and ATc are proportional to | J-|^, the square of 
the negative exchange parameter leading to the Kondo 
effect. Dy’s volume collapse itself has been suggested to 
have its origin in the Kondo volume collapse model of 
Allen and Martin .1^ 

Does perhaps the same scenario apply for Tb? In Fig. 
6 the pressure dependence of the superconducting transi¬ 
tion temperature To{P) of the dilute magnetic alloy Y(0.5 
at.% Tb) from our previous worlP is compared to that 
for elemental Y metal.l^ To a pressure of ~50 GPa, To 
for the dilute magnetic alloy is seen to increase with pres¬ 
sure at the same rate as for Y. However, just above the 
pressure of Tb’s volume collapse at 53 GPa, the Tc{P) 
dependence for the alloy begins to pull away rapidly from 
that of Y, reaching a maximum suppression ATc ~ 5 K 
at 81 GPa, the highest pressure of the experiment. This 
strong suppression of Y’s superconductivity by dilute Tb 


TABLE I. Values for Tb of the average To and maximum 
magnetic ordering temperatures, spin-disorder resis¬ 
tance for T > and phonon factor a as a function 

of pressure from runs 1 and 2 (see text). 


run 

P(GPa) 

rc(K) 

ymax(K) 

RTrimn) 

a 

1 

2 

211 

214 

881 

1.16 

1 

5 

117 

163 

523 

1.24 

1 

8 

73 

84 

380 

1.24 

1 

18 

64 

70 

154 

1.0 

1 

27 

95 

145 

513 

0.66 

1 

34 

101 

165 

721 

0.72 

1 

44 

112 

188 

769 

1.02 

1 

53 

105 

196 

746 

1.11 

1 

61 

101 

179 

518 

1.05 

1 

68 

86 

154 

479 

0.89 

1 

77 

72 

105 

469 

0.74 

1 

88 

73 

145 

1001 

0.69 

1 

97 

83 

173 

1140 

0.69 

1 

109 

107 

206 

1165 

0.69 

1 

114 

131 

230 

1189 

0.69 

1 

119 

150 

253 

1201 

0.69 

1 

126 

172 

291 

1211 

0.69 

1 

135 

192 

323 

1214 

0.69 

2 

2 

208 

210 

83 

0.98 

2 

5 

103 

142 

43 

0.98 

2 

9 

57 

66 

36 

0.98 

2 

16 

52 

66 

13 

1.0 

2 

35 

106 

177 

111 

0.58 

2 

52 

111 

199 

138 

0.92 

2 

61 

101 

183 

103 

0.79 

2 

73 

82 

159 

91 

0.56 

2 

86 

69 

113 

152 

0.41 

2 

94 

76 

154 

219 

0.41 

2 

99 

88 

180 

234 

0.41 

2 

106 

111 

205 

239 

0.41 

2 

112 

133 

232 

247 

0.41 

2 

116 

148 

249 

250 

0.41 

2 

120 

163 

274 

252 

0.41 

2 

126 

183 

306 

252 

0.41 

2 

141 

203 

341 

243 

0.41 


ions points to giant Kondo pair breaking, as has previ¬ 
ously been observed in high pressure studies o n th e di¬ 
lute magnetic alloys L^Ce),I^La(Pr),l^ Y(Pr),®^and, 
most recently, Y(Dy).® In contrast, as seen in the in¬ 
set to Fig. 6, Tc{P) for Y(0.5 at.% Gd) does not begin 
to deviate markedly from that of Y metal near 59 GPa, 
where Gd’s volume collapse occurs, but rather faithfully 
tracks Y’s value of T^ to 127 GPa, the maximum pres¬ 
sure of the experiment. Unlike for Tb, the magnetic state 
for Gd ions in Y remains stable to 127 GPa, so that no 
Kondo phenomena are expected. We thus suggest that 
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FIG. 6. (color online) Tc versus pressure for Y(0.5 at.% Tb) 
compared to that for Y, inset showing similar graph for Y(0.5 
at.% Gd)Vertical dashed line marks pressure of volum e col- 
lapse for Tb at 53 GPrP and in inset for Gd at 59 GPa.^^^^ 
At top of graph are crystal structures taken on by supercon¬ 
ducting host Y.l^ 

the anomalous pressure dependences To{P) and ATc{P) 
in Tb and Y(Tb) alloy, respectively, have their origin in 
Kondo physics, as does Tb’s volume collapse itself. In 
support of these suggestions we point out that XANES 
and XES experiments on Tb to extreme pressure reveal 
that neither a change in valence nor a magnetic local- 
itinerant transition occur to a pressure of ~ 70 GPa, well 
above the volume collapse pressure for Tb at 53 GPa.l^ 
Could perhaps an alternative explanation for the 
anomalously high magnetic ordering temperatures To in 
Tb be the effect of crystalline electric fields? It has 
been argued that such fields are likely responsible for 
the significant enhancement of To over de Gennes scal¬ 
ing in a s eries of i?Rh 4 B 4 compounds, where i? is a 
lanthanide If th e magnetic anisotropy is strong, it 
has been showipil^SI that the crystal field enhancement 
can be as large as the factor 3Jt/{Jt +1) = 2.6 for triva- 
lent Tb where T = 3, S' = 3, and Jt = 6. No crystal 
field effects are possible for Gd since it carries no or¬ 
bital moment (T = 0). The lack of a sharp upturn or 


other anomalies in To and ATc in the pressure region 60 
- 127 GPa would be consistent with the certain absence 
of crystal field effects in Gd. The fact that the pres¬ 
sure dependence of To is very similar for both Gd and 
Tb to 53 GPa indicates that crystal field effects in Tb, if 
present, are only significant for pressures above 53 GPa 
where the To(T) dependence becomes anomalous. In a 
crystal field scenario, however, it would be difficult to 
understand the sharp upturn in the suppression of su¬ 
perconductivity ATc in the dilute magnetic alloy Y(0.5 
at.% Tb) for pressures above 53 GPa. This strong sup¬ 
pression of superconductivity points rather to a Kondo 
physics scenario with strong Kondo pair breaking. 

Further experimentation is necessary to unequivocally 
establish the origin of the anomalous behavior of To and 
ATc in Tb and Y(Tb) alloy, respectively, for the pres¬ 
sure region above 53 GPa. Such experiments could in¬ 
clude an extension of the pressure range to 2 Mbar to 
search for the characteristic “Ko ndo sinkhole behavior” 
in Tc(T) observed for Y(Pr),l^ La(Ce),l^ and La(Pr)P 
where the Tc-suppression ATc reaches a maximum as the 
Kondo temperature Tk passes through the experimental 
temperature range, but falls off again at higher pressures 
where Tk far exceeds T^. Inelastic neutron or x-ray scat¬ 
tering studies to extreme pressures would help establish 
whether crystal-field splittings play a role in the anoma¬ 
lously high values of To for Tb. 

In summary, measurements of the electrical resistivity 
of Tb metal to extreme pressures reveal that the magnetic 
ordering temperature To exhibits a highly non-monotonic 
pressure dependence, appearing to rise for P > 80 GPa 
to anomalously high values. Parallel experiments on Gd 
and dilute magnetic alloys of Gd and Tb with Y suggest 
that under extreme pressures Tb is transformed from a 
magnetically conventional lanthanide into one with an 
unconventional magnetic state, perhaps involving Kondo 
physics, with anomalously high values of Tq. In contrast, 
Gd remains a magnetically conventional lanthanide to 
pressures of at least 127 GPa. 

Acknowledgments. The authors would like to thank 
T. Matsuoka and K. Shimizu for sharing information on 
their high-pressure electrical resistivity techniques used 
in the present study. This work was supported by the 
National Science Foundation (NSF) through Grant No. 
DMR-II04742 and by the Garnegie/DOE Alliance Cen¬ 
ter (CDAC) through NNSA/DOE Grant No. DE-FC52- 
08NA28554. Work at Argonne National Laboratory is 
supported by the US Department of Energy, Office of 
Science, under contract No. DE-AC02-06CHI1357. 


Current address: Department of Condensed Matter 
Physics and Material Science, Brookhaven National Labo¬ 
ratory, Upton, NY 11973, USA 
jss@wuphys.wustl.edu 


J. Lim, G. Fabbris, D. Haskel, and J. S. Schilling, Phys 
Rev. B 91 , 045116 (2015). 

See: K. N. R. Taylor and M. I. Darby, Physics of Rare 
Earth Solids (Chapman and Hall Ltd., London, 1972). 


















® W. C. Koehler, J. Appl. Phys. 36, 1078 (1965). 

^ D. B. McWhan and A. L. Stevens, Phys. Rev. 139, A682 
(1965). 

® D. D. Jackson, V. Malba, S. T. Weir, P. A. Baker, and Y. 

K. Vohra, Phys. Rev. B 71, 184416 (2005). 

® M. Mito, K. Matsumoto, Y. Komorida, H. Deguchi, S. 
Takagi, T. Tajiri, T. Iwamoto, T. Kawae, M. Tokita, K. 
Takeda, J. Phys. Chem. Solids 70, 1290 (2009). 

^ S. A. Thomas, J. M. Montgomery, G M. Tsoi, Y. K. Vohra, 
G. N. Chesnut, S. T. Weir, C. A. Tulk, A. M. dos Santos, 
High Press. Res. 33, 555 (2013). 

® G. Fabbris, T. Matsuoka, J. Lim, J. R. L. Mardegan, K. 
Shimizu, D. Haskel, and J. S. Schilling, Phys. Rev. B 88, 
245103 (2013). 

® N. C. Gunningham, W. Qiu, K. M. Hope, H.-P. Liermann, 
and Y. K. Vohra, Phys. Rev. B 76, 212101 (2007). 
Material Preparation Genter, Ames Lab, US DOE, Ames, 
Iowa http:/www. mpc.ameslab.gov. 

James S. Schilling, Mat. Res. Soc. Symp. Proc. 22, 79 
(1984). 

K. Shimizu, K. Amaya, and N. Suzuki, J. Phys. Soc. Jpn. 
74, 1345 (2005). 

W. B. Daniels and W. Ryschkewitsch, Rev. Sci. Instr. 54, 
115 (1983). 

A. D. Ghijioke, W. J. Nellis, A. Soldatov, and I. F. Silvera, 
J. Appl. Phys. 98, 114905 (2005). 

Y. Akahama and H. Kawamura, J. Appl. Phys. 100, 
043516 (2006). 

S. Klotz, J. S. Schilling, and P. Muller, in Frontiers of High 
Pressure Research, edited by H. D. Hochheimer and R. D. 
Etters (Plenum, New York, 1991) p. 473. 

M. Debessai, J. J. Hamlin, and J. S. Schilling, Phys. Rev. 
B 78, 064519 (2008). 


B. R. Coles, Adv. Phys. 7, 40 (1958). 

R. V. Colvin, S. Legvold, and F. H. Spedding, Phys. Rev. 
120, 741 (1960). 

R. Patterson, C. K. Saw, and J. Akella, J. Appl. Phys. 95, 
5443 (2004). 

H. Hua, V. K. Vohra, J. Akella, S. T. Weir, R. Ahuja, 
and B. Johansson, Rev. High Pressure Sci. Technol. 7, 233 
(1998). 

22 D. Errandonea, R. Boehler, B. Schwager, and M. Mezouar, 
Phys. Rev. B 75, 014103 (2007). 

2® J. Duthie and D. Pettifor, Phys. Rev. Lett. 38, 564 (1977). 
21 G. S. Fleming and S. H. Liu, Phys. Rev. B 2, 164 (1970); 

S. H. Liu, Phys. Rev. 127, 1889 (1962). 

2® J. W. Allen and R. M. Martin, Phys. Rev. Lett. 49, 1106 
(1982). 

2® Z. P. Yin and W. E. Pickett, Phys. Re. B 74, 205106 
(2006). 

21” B. T. Matthias, H. Suhl, and E. Corenzwit, Phys. Rev. 
Lett. 1, 92 (1958). 

2® M. B. Maple, Appl. Phys. 9, 179 (1976). 

2® G. K. Samudrala, G. M. Tsoi, and Y. K. Vohra, J. Phys.: 
Condens. Matter 24, 362201 (2012). 

J. Hamlin, V. Tissen, and J. S. Schilling, Physica C 451, 
82 (2007). 

M. Maple, J. Wittig, and K. Kim, Phys. Rev. Lett. 23, 
1375 (1969). 

®2 J. Wittig, Phys. Rev. Lett. 46, 1431 (1981). 

J. Wittig, Valencies Instabilities, edited by P. Wachter and 
H. Boppart (North-Holland, Amsterdam, 1982) p. 427. 

D. R. Noakes and G. K. Shenoy, Phys. Lett. 91A, 35 
(1982). 

B. D. Dunlap, L. N. Hall, F. Behroozi, G. W. Crabtree, D. 
G. Niarchos, Phys. Rev. B 29, 6244 (1984). 



